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5 ■ Abstract 

-p^ ' A calculus for the derivatives of the eigen values of the neutralino mass matrix 

Ph| with respect to the CP violating background fields is developed and used to com- 

Q-i. pute the mixings among the CP even and the CP odd Higgs sectors arising from 

^ . the inclusion of the neutralino sector consisting of the neutralino, the Z boson, and 

' the neutral Higgs bosons {y^ — Z — h^ — H^) exchange in the loop contribution to the 

^ ' effective potential including the effects of large CP violating phases. Along with 

5t ! the stop, sbottom, stau and chargino-W-charged Higgs {x^ — W — H^) contribu- 

tions computed previously the present analysis completes the one loop corrections 
to the Higgs boson mass matrix in the presence of large phases. CP violation in the 
neutral Higgs sector is discussed in the above framework with specific focus on the 
mixings of the CP even and the CP odd sectors arising from the neutralino sector. 
It is shown that numerically the effects of the neutralino exchange contribution on 
the mixings of the CP even and the CP odd sectors are comparable to the effects 
of the stop and of the chargino exchange contributions and thus the neutralino 
exchange contribution must be included for a realistic analysis of mixings in the 
CP even and the CP odd sectors. Phenomenological implications of these results 
are discussed. 



1 Introduction 



CP violation in supersymmetric theories via soft supersymmetry breaking param- 
eters has received a considerable degree of attention since the beginning of the 
formulation of supersymmetric models]^]. Recently, there has been enhanced in- 
terest in the investigation of their effects due to the realization that supersym- 
metric theories may allow for large CP violating phases]^ consistent with the 
electric dipole moment of the electron and of the neutron^]. Such a situation can 
arise because of several possibilities, such as the SUSY spectrum being heavy 0], 
due to internal cancellations [|] and due to the possibihty that the CP phases 
may reside in the third generation and consequently their effects on the first two 
generation EMDs are suppressed 0. Of course it is possible that a more uni- 
fied framework may determine the combination of phases that enter the EDMs 
to be small|p. However, we shall investigate here the possibility that the phases 
are large and the EDM constraints are satisfied by one of the methods discussed 
above so that the sparticle spectrum is consistent with the naturalness constraints 
(see, e.g., Ref. 0). In this case their effects on low energy physics can be quite 
significant and a number of low energy phenomena have been discussed includ- 
ing the effect of CP phases. These include the effect of CP phases on g-2P], on 
dark matter on the trileptonic signal on baroyogenesis [p!l|] , and on other 
low energy phenomena||12||. Another area where the effect of CP phases has been 
discussed is the Higgs sector |jl3|, |14|, [T^, |16|, II^, |l8i. It is well known that loop cor- 



rections to the Higgs masses and mixings are important ||T9[. In fact in the absence 
of the loop corrections the lightest Higgs boson mass must lie below Mz which is 
already exprimentally excluded and it is the presence of the loop corrections that 
raises its value above Mz- An interesting phenomenon arises if the loop corrections 
have CP violating phases. In this case it has been pointed out that a significant 
mixing can occur between the CP even and the CP odd neutral Higgs sectors of 
the theory 0. In Refs.[|l|, 0, |16|, [l^, [T|] the effect of CP phases via the stop 
and sbottom exchanges was carried out. Further, in the work of Ref. |T6[ it was 
pointed out that the effect of chargino loop corrections can be quite significant 
and in fact the CP effects from the chargino exchange may even dominate the CP 
effects from the stop-sbottom exchange for the case of large tan (3. 

In this paper we give an analysis of the one loop correction to the Higgs boson 
mass including the neutralino-Z boson-neutral Higgs exchange including the CP 
violating phases. The inclusion of the CP dependent neutralino exchange correc- 
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tions are more intricate relative to the stop-sbottom exchanges and the chargino 
exchanges. This is due to the fact that the stop-sbottom exchange and the chargino 
exchange involve diagonalization of only 2x2 squark and chargino mass matrices 
and thus the evaluation of their contribution can be carried out analytically in a 
straightforward fashion. For the case of the neutralino exchange the neutralino 
mass matrix is a 4 x 4 object and its diagonalization analytically is more intricate 
and a straightforward technique for the analysis is wieldy. In this paper we develop 
a calculus for the derivatives of the eigen values of the neutralino mass matrix to 
obtain an explicit analytic expression for the neutralino exchange contribution. 
The outline of the rest of the paper is as follows. In Sec. 2 we give the Higgs poten- 
tial and discuss the minimization conditions in the presence of the CP violating 
phases. In Sec. 3 we discuss the calculus for the computation of derivatives of the 
eigen values of the neutralino mass matrix. In Sec. 4 we use the technique of Sec. 3 
and compute the one loop contributions to the Higgs boson mass matrix from the 
neutralino- Z-neutral Higgs boson exchange. Discussion of the numerical results is 
given in Sec. 5. Conclusions are given in Sec. 6. Some further details of the analysis 
are given in Appendices A and B. 



2 CP Phases and Minimization of Higgs Poten- 
tial 

We begin by defining the soft SUSY breaking parameters for the mSUGRA casepO|. 



Here the low energy physics for the CP conserving case is parametrized by mg, 
mi, Aq, and tan/3 where mo is the universal scalar mass, mi is the universal 

2 2 

gaugino mass, Aq is the universal trilinear coupling, and tan (3 = ^ is the ratio 
of the Higgs VEVs, where the VEV of H2 gives mass to the up quarks and the 
VEV of Hi gives mass to the down quarks and the leptons. In the presence of CP 
violation mSUGRA allows for only two CP violating phases which can be taken 
to be and a^o where is the phase of the Higgs mixing parameter /iq and 
is the phase of Aq. The analysis of this paper, however, will be more general, 
valid for the MSSM parameter space. The Higgs sector in MSSM at the one loop 
level is described by the scalar potential V{Hi, H2) = Vq + AV where 



Vq = ml\Hi\^ + ml\H2\^ + {mlHi.H2 + H.C.) 

— \H H P I ~ 
2 I 1- 21 ^ 



l^i±^lH- |4 , M±^lH- |4 9I rr rr \2 .ial-ai), rr \2\rT \2 
7^ "I 7^ l-n2| ;:rl-"l--"2| H |-ni| |-n2| 
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Here ml = m?^^ + \^\'^^m\ = m?^^ + |/ip,m| = \^B\ and mu^^ and -B are the 
soft SUSY breaking parameters, is the one loop correction to the effective 
potential pi], ^ and includes contributions from all the fields that enter MSSM 



consisting of the standard model fields and their superpartners, i.e., the sfermions. 



the gauginos and higgsinosp2|. The sum over i in Eq.(l) runs over particles with 
spin Ji and Ci{2Ji + 1) counts the degrees of the ith particle, and Q is the renor- 
malization group running scale. It is well known that the one loop corrections 
to the effective potential can make significant contributions to the Higgs vacuum 
expectation values in the minimization of the effective potential . 



In general the effective potential depends on the CP violating phases and its 
minimization will lead to induced CP violating effects on the Higgs vacuum expec- 



tation values[0]. It is found convenient to parameterize the Higgs VEVs in the 
presence of CP violating effects in the following form 

/ m 



where Oh is in general non-vanishing as a consequence of the minimization condi- 
tions. Thus the minimization of the potential with respect to the fields 0i, ^1, 02, ^2 
gives 



1 2 . a 

— (— — )o = ^3 smOn 

^ )o = + (f 1 — f 2 ) + tan [3 cos 6^ (3) 



and 



1 .^^^N 2 . a 

— (— — )o = rngSmfc^H 

Vi d%l)2 

(^:^)o = "^2 z — [v^-v^) + m^cotp cos Oh (4) 

V2 0<p2 8 

In the above the subscript stands for the fact that we are evaluating the relevant 
quantities at the point 0i = 02 = "^i = "02 = 0. We note in passing that in Eqs.(3) 
and (4) only one of the two equations that involve the variation with respect to 
ipi and '02 is independent 
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3 Calculus for Derivatives of Eigen Values of Neu- 
tralino Mass Matrix 



As mentioned in Sec.l, in previous analyzes computations of the CP dependent 
loop corrections from the stop-sbottom and from the chargino- W- charged Higgs 
sectors have been carried out. In these analyzes one was able to analytically obtain 
the eigen values by diagonalizing the 2x2 squark matrices and the 2x2 charginio 
mass matrix and then differentiate them analytically to obtain the loop correction 
to the Higgs mass matrix. As also pointed out in Sec.l for the neutralino exchange 
case the situation is more difficult since the neutralino mass matrix is a 4 x 4 matrix 
and the analytic solutions for the eigen values of the neutralino (mass)^ matrix 
are not easily obtained. Here we expand on a technique introduced in Ref. to 
derive a calculus for the derivatives of the eigen values for the neutralino mass 
matrix. This technique is valid for an arbitrary high order eigen value equation. 
We shall show that quite remarkably even though one cannot analytically solve 
for the eigen values one can analytically solve for the derivatives of the eigen 
values with respect to the background fields in terms of the eigen values and the 
parameters that appear in the eigenvalue equation. To illustrate the procedure we 
consider an nth order eigen value equation 



F(A) = Det{M^M - XI) = A" + c^^-^U""^ + c("-2)A"-2 + .. + c^^^A + c^") = (5) 
Here the co-efficients are explicit functions of the background fields 

= {01, 02,^1, V'2} (6) 

while the eigen values are implicit functions of the background fields through the 
satisfaction of the eigen value equation. Eq. (5) has n eigen values which we denote 
by Aj {i = 1,2, ..,n). From Eq.(5) it follows that 

dXi ^DaF 

)a=a, [ij 



and 



d^Xi , D^FDsFDlF , D^FDsD^F + DpFD^DxF D^D.F 



H TW-P^ A=A, 



d^^d^p ' (DxF)^ (DxF)^ DxF 

where D\ differentiates the A dependence in F 

D.FW = § (9) 



4 



and Da differentiates only the co-efficients in Eq.(5), i.e., 

DaF = ci"-i)A("-^) + ci"-2)A("-2) + .. + c«A + ci°) (10) 
DoiDj^F are similarly defined where c^'^^ etc are replaced with c^p where 

and the derivatives DaD\ are defined in an obvious way. We note in passing that 
Da and Dx commute 

[Da,Dx] = (12) 

Eqs.(7) and (8) are the central equations of our analysis. It is easy to check that 
for the 2x2 matrix case, e.g., for the stop and the chargino exhanges, they give 
exactly the results gotten by explicit differentiation of the eigen values. However, 
now these equations provide us with a technique of analyzing cases where the 
analytic solutions to the eigen values are not available. 

4 Neutralino, Z and neutral Higgs loop contri- 
butions 

As mentioned above the CP dependent contributions to the Higgs boson masses 



from stop and sbottom exchanges have been discussed at length in the literature |13 



Tj, |T5|, |T6|, |l^, [T^. More recently the CP dependent chargino- W-charged Higgs 



contributions were also discussed [T^]. In this work we use the technique discussed 



in Sec. 3 to compute the contribution from the neutralino -Z - neutral Higgs ex- 
change. The loop correction in this sub sector is given by 



647r2 



4 M^o M2 



+M^olog{ 



3. 
2' 



- -) + 



HO 



(13) 



The neutralino mass matrix is given by 



/ 


nil 







V2^2 


\ 







1712 


92 ttO 


92 ttO 
72^2 













-/^ 




I 


9i_ZjO 
v^^2 


92 ttO 
72^2 


-/^ 





/ 



(14) 
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where /i = \fi\e^^'^, rhi = |mi|e*^^ and m2 = |m2|e*^^. We note that in the super- 
symmetric hmit M^o = (0, 0, Mz, Mz) and (M/jO, Mhq) = {Mz, 0) and consequently 
in this hmit the loop corrections from this sub sector vanish. We return now to the 
full analysis and follow the method described in Ref . [|T6| to minimize the potential 
and compute the loop corrections. First we give the determination of 6h from the 
minimization constraints including the stop, the sbottom, the stau, the chargino 
and neutralino contributions. One finds that 6h is given by the equation 



mlsm9H = ^/3/,J/i||A| sin7i/i(m|,m|) + ^/5?,J/i| |A| sin 7^/1 (m?^,m?J 



1 ^2 
+ ^(^h^f^\\Ar\ sin7^/i(m?^,m?J - J^lf^W^^l sin72/i(m|^, m|^, 

4 

up 

+ (— (72l^^inyu|^|^2| sin 72 — (y'^|m2p|/i|^|mi| sin 71)) (15) 



I 4 M^o 

H ^T^(^^(~7#) ~ l)(^x°(~^2l/^ll'^2| sin 72 - gl\iJ.\\mi\ sin 71) 



+Mj,(5'2(l"^iP + lAini"^2||/i| sin 72 + fi-idmip + |/ip)|mi||/i| sin 71) 



where 



Dj = (DaF)a=a, = 4Mjo + SaM^o + 2bM^o + c 

Pht — TTT^J Pht — TTT^iPhr 



167r2' 167r2"^"^ 16772 

7t = oiAt + lb = oiAi + 6*/., 7r = oiAr + d^,, 7i = ^1 + 6^^, 72 = 6 + (^6) 

and where a,b,c are defined in Appendix A and fi{u,v) is given by 

/i(n, t;) = -2 + /o(?— + log- (17) 

(^^ f — M U 

To construct the mass squared matrix of the Higgs scalars we need to compute the 
quantity 



where M^^p is the contribution from Vq and AM^^ is the contribution from AV 
where $0(0; = 1 — 4) are defined by Eq.(6) and as already mentioned earlier the 
subscript means that we set 0i = 02 = V^i = ^2 = after evaluating the mass 
matrix. The loop contribution AM^^ arising from the neutralino- Z-neutral Higgs 
sector is given by 
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AAe, = ^.Str(^g^los^ + M=^^(,„,_ - 1))„ (19) 
Computation of the 4x4 Higgs mass matrix in the basis of Eq.(6) gives 



-{Ml + Ml) s pep + Ai2 MI4 + MI4 + A22 A23S/3 A23CP 

AisSp A23SP (Ml + A33)4 (Ml + A33).W 

V Ai3C^ A23C^ (Ml + A33)S/3C/3 (Ml + A33) 4 / 

(20) 

where 0/3(5^3) = cos/3(sin/3) and is given by 

m\ = (sin/3cos/3)~^(-m3COS^ + ^/3ftJA||/^|cos7t/i(m? ,m|) 

1 1 

+ 2(^h,\Ab\\n\cos ^bfi (m?^ , + 2 /^/^r I 1 1 I cos 7^/1 (m?^ , m?^ ) 

+Tfel^2||/i|cos72/i(m^,m\) 




+M^ (5(^(|mi|2 + l/ini/illmal cos 72 + gl{\m2\^ + |/ip) |/i| |mi| COS71) 



— 5f2l^inA*l^l^2| COS 72 — g'i|m2nAi|^|?7ii| COS71]) (21) 

The first term in the second brace on the right hand side of Eq. (21) is the tree term, 
while the second, the third, the fourth and the fifth terms come from the stop, 
sbottom, stau and chargino exchange contributions. The remaining contributions 
in Eq.(21) arise from the neutrahno sector. The A's appearing in Eq.(20) can be 
decomposed as follows 

Aa/3 = A^0 + A^pi + A^pf + A^i3^+ + A^p^o (22) 

where A^^i is the contribution from the stop (and top) exchange in the loops, A^^^ 
is the contribution from the sbottom (and bottom) exchange in the loops, A^pf 
is the contribution from the stau (and tau) exchange, Aap^+ is the contribution 
from the chargino (and W and charged Higgs) exchange in the loops, and A^p^A) is 
the contribution arising from the neutralino (and Z and neutral Higgs exchange) 
in the loops. The computations of A^^i, A^^j^, A^isf, and A(^^^+ have been given 
before and would not be reproduced here. We compute here only the Aq^^o arising 
from the (x° — Z — — H^) exchange. The A^^^o are listed below. 
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(aiM^. + biM^. + ciMl + di)\l2M^. + daM^ + 2b) 
^2(aiM6. + biM^. + ciMl + di){3aiM^. + ^b^M^. + ci) ^ 



where £)j is defined in Eq.(16) and /2 is defined by 

f,{u,v)^-2 + ^ln- (24) 

V — M M 



(asM^. + fcsMj . + C2M2. + d2f{l2M^. + GaM^. + 26) 
^2{a2Ml. + + C2M2. + d2)(3a2M^. + 262M2. + c^) ^ 
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2^ ^2i?0 M^o 



1 4 j^^2 

A,..o = -^EMi(M^)-i) 

(aiM^. + 6iM^. + ciM2. + rfi)(a2M6. + ^aM^. + csM^. + rf2)(12M^. + GaM^. + 26) 
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(aiM^. + biM^. + c,M^. + rfi)(3a2M^. + 2b2M^. + C2) 

(a2M6. + + C2M2. + d2){3a,M^. + 2b,Ml + ci) 

1 ^ {a,Ml + + c^Ml + d,){a2Ml + 62M4. + C2MI + 6^2) , il^ 

IGtt^ g2 ) 

3 



■(77 r.. 2 ..2^9 /2«o, M,^o) - -(^f + ^2') V^2/n- 



4 



1 / 2 , 2n t'l5o + 'y2^0 , ^'^I'vofi^ 

-Iefa+g^) (Mi.-M|.) '"M|:)P«' 

Ai3xO = -T7-iE<.(MTf)-l)- 



167r%tl <9 sin/? 

(aiM^. + biM^. + CiM2. + di)(a3M6. + b^M^. + csM^ + d3){12M^. + eaM^ + 2b) 

^ (aiM^. + b,M^. + CiM^. + c^i)(3a3Mj. + 263M^. + cg) 



^1 

+ ^2 I 



J_ f 1 (aiM^^. +b,M^^^. +c,M^^. +dO(a3M^^^. +b3M^^^. +C3M^^. +d3) M 

(a;,M6. + 6;,M^. + 4M^. + 4)(a2M^. + b2M^^ + C2M^. + d2)(12M^4^. + GaM^ + 26) 

{a',Ml + 6;,M^4^ + 4M^^. + d',)i3a2M^^ + 262M^^ + C2) 
(a2M^6. + 62M^4. + C2MI + ci2)(3a-3M^4^. + 2b',Ml + 4) 

167r%tlcos/3 D| '^lg2A^«i 

A33.o = -^|:Mi.(M^)-l) 
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( 



{a',Ml + h',M^^ + d^Ml + 4)'(12M^. + 6aM^. + 26) 1 

D| cos2 13 

2{a',Ml + 6- M^^ + c'.Ml + 4) (Sa^M^^^. + 2b',Ml + 4) 1 

cos2 /3 



= —]n(^^) (2()) 

167r%tlcos2/3 D] ^""^ ' 

The parameters a,b,c and the derivatives ai,bi,Ci,di (i=l,2, etc.) that appear in 
Eqs. (23-29) are defined in Appendices A and B. Eqs. (23-29) constitute the main 
new theoretical results of this paper. These results along with the computations of 
A^^g, Aai3f and Aai3x+ ^ complete determination of the CP dependent 
one loop contributions to the Higgs boson masses and mixings. As has been noted 
before it is preferable to work with a 3 x 3 matrix rather than the 4x4 matrix 
of Eq.(20). The desired 3x3 matrix can be gotten from Eq.(20) by going to the 
basis 

ijjiD = sin /St/ji + cos /?-?/'2, i>2D = - COS P^pi + sin /3i'2 (30) 

In this basis the field ■02£i is the zero mass Goldstone boson and decouples while 
the remaining {massy matrix in the basis 0i, 02, i^iD is given by 

M|4 + Ml4 + Aii -(M| + Ml)s^c^ + Ai2 A13 \ 
-(M| + Ml)spcp + A12 M|4 + + A22 A23 

Ai3 A23 {Ml + A33) / 

(31) 

We label the eigen values for this case "m^^j ""^^3 corresponding to the eigen 

states Hi, H2, H-^. These eigen states are in general admixtures of the CP even 
and the CP odd states due to the mixing generated by A13 and A23. Thus the CP 
even-odd mixings arise from A13 and A23 and these are nonvanishing only in the 
presence of CP violation and vanish when the phases go to zero and one recovers 
the usual result of two distinct (one CP even and the other CP odd) Higgs sectors. 
We note in passing that A33 also vanishes in the limit when the CP phases go 
to zero. This was also the behavior that was observed when the contributions 
from the stop, sbottom, stau and chargino exchanges were considered. Since the 
main point of this work is to study the phenomenon of CP even-odd mixing the 
main focus of our analysis is the computation of Aij and specifically of A13 and 
A23 which are the basic sources of mixings between the CP even and the CP odd 
sectors. We order the eigen values of Eq.(31) in such a way that in the limit 
of no CP violation one has (m//^, mj^j, ""^Hs)— {^H,^h,^A) and {Hi, H2, H^)^ 
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{H, h, A) where (h, H) are (light, heavy) CP even Higgs and A is the CP odd Higgs 
in the absence of CP violation. 

5 Discussion of the Neutralino Exchange Con- 
tribution to CP even CP odd Higgs Mixing 

The analytical results given above are quite general as they apply to the MSSM 
parameter space. However, the MSSM parameter space is quite large. Thus for 
a numerical study of the CP effects including those from the neutralino sector we 
will work with a constrained set of parameters consisting of the parameter space 
mo, mi, niA, \Ao\, tan/?, aAo, ^i, ^2 and ^3. Starting with these all other low 
energy parameters are obtained by a renormalization group evolution by running 
the parameters from the GUT scale down to the electro-weak scale. Of course one 
is free to utilize the formulae derived above for the more general MSSM parameter 
space. As discussed in Sec.l one can satisfy the EDM constraints in the presence 
of large phases. This can come about in a variety of ways. As pointed out in Sec.l 
one possibility is that the internal cancellations can occur which allow for large 
phases consistent with the EDM constraints. The other possibility is that that CP 
phases appear only in the third generation which suppresses their contributions 
to the EDMs of the quarks and the leptons in the first two generations to achieve 
consistency with the experimental constraints. There also exist scenarios which are 
linear combinations of these two. For the purpose of this analysis we do not revisit 
the problem of the satisfaction of the EDM constraints. Rather we shall assume 
that regions of the parameter space exist where such constraints are satisfied and 
examine the effect of the phases on the Higgs masses and mixing. Specifically we 
are interested in the effects of the neutralino exchange contributions on A13 and 
A23, and thus, on the mixings of the CP even and the CP sectors. 

It was pointed out in Sec. 4 that the neutralino, the Z and the neutral Higgs 
exchanges together form a sub sector so that in the supersymmetric limit one 
finds that the one loop correction to the effective potential from this sub sector 
vanishes. This phenomenon is similar to what was also seen in the exchange of the 
chargino, the W and the charged Higgs where the contribution from that sector to 
the effective potential vanishes in the supersymmetric limit. It was also seen in the 
analysis of the chargino- W-chargcd Higgs exchange that the CP even-odd mixing 
arising from this sector was roughly Q independent because of the sum of the three 
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separate contributions within this sector. A very similar situation is also realized in 
the neutralino sector. Here again because of the contributions from the neutralino, 
the Z and the neutral Higgs exchanges their sum contribution to the CP even-odd 
mixing is roughly scale independent. However, unlike the chargino-W- charged 
Higgs exchange where one could demonstrate the above phenomenon analytically, 
here one has to demonstrate it numerically due to the more analytically complex 
nature of the results. This is exhibited in Fig.l where a plot the percentage of the 
CP even component 0i and the CP odd component ipiD of Hi as a function of Q is 
given. The analysis shows an approximate independence in Q of the CP even-odd 
mixing. We turn now to a discussion of other aspects of the analysis below. 

In Fig. 2 we plot the quantity A13 as a function of the CP phase of the U(l) 
gaugino mass .^1. The plots exhibited in Fig.2 contain the stop, the sbottom, the 
stau the chargino and the neutralino exchange contributions. Among the above 
exchanges the neutralino exchange contribution is the only one that depends on ,^1, 
and thus the variation of A13 with ^1 arises only from this exchange. From Fig.2 
the size of the neutralino exchange contribution can be seen to be fairly substantial. 
Specifically, the analysis of Fig.2 shows that the neutralino exchange contribution 
is comparable to the effects from the stop and chargino exchanges. A plot of A23 
vs ^1 is given Fig.3. As in Fig.2 one finds that A23 is quite sensitive to the CP 
violating phase ^1. As in Fig.2 here again the neutralino exchange contribution 
is comparable to the stop and the chargino exchange contribution. An analysis 
of the percentage of the CP even component 0i of Hi (upper curves) and of the 
percentage of the CP odd component ipiu of Hi (lower curves) arising from the 
exchange of the stop, the sbottom, the stau, the chargino and the neutralino sector 
contributions as a function of ,^1 is given in Figs. 4. As expected from the analysis 
of Fig.2 and Fig.3 one finds that there is a significant mixing between the CP even 
and the CP odd components of Hi. Further, as also expected from the analysis 
of Figs 2 and 3, the CP even and CP odd components of Hi show a reasonably 
strong dependence on ^1. 

An analysis of the CP even and CP odd mixing in Hi as a function of the SU(2) 
gaugino phase is given in Fig.5. Unlike Figs. 2-4, where the entire ,^1 dependence 
arose from the neutralino exchange contribution here the .^2 dependence of the CP 
even and CP odd components of Hi arises from two sources, i.e., from the chargino 
and the neutralino exchange contributions. Because of this the dependence of the 
CP even and CP odd components on {2 is much stronger than on ^1 as may be 
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seen by comparing the plots of Figs. 2-4 with the plots of Fig. 5. In Fig. 6 a plot 
of the percentage of the CP even component 0i of Hi (upper sets) and the CP 
odd component ipiD of Hi (lower sets) arising from the exchange of the stop, the 
sbottom, the stau, the chargino and the neutralino sector contributions is given 
as a function of 6'^. In this case we find that the dependence of the CP even and 
the CP odd components on 9^ is also very strong. Indeed in this case the mixings 
between the CP even and the CP odd states can be maximal depending on the 
value of 0^. The strong dependence on 0^ can be understood as due to the fact 
that all contributions, i.e., the stop, the sbottom, the stau, the chargino, and the 
neutralino contributions, depend on This in contrast to the dependence on 
which arises only from the neutralino exchange. 

Finally, in Fig. 7 we give an analysis of the percentage of the CP even com- 
ponent 01 of Hi (upper sets) and the CP odd component ipin of Hi (lower sets) 
arising from the exchange of the stop, the sbottom, the stau, the chargino and the 
neutralino sector contributions as a function of tan (3. We find that the CP even 
and the CP odd mixings show a strong dependence on tan (3. A similar strong 
dependence on tan/5 was seen also in previous analyses|T^. We note that the in- 
clusion of the neutralino contribution further sharpens the tan (3 dependence and 
one finds that the CP even (odd) component can vary from 100% (0%) to less 
than 60% (more than 40%) as tan (3 is varied. This sharper behavior of the am- 
plitudes with tan/5 arises from the additional contributions from the neutralino, 
the neutral Higgs and the Z boson exchanges. An analysis similar to the above 
can be carried out for the case of the H2 and H^ fields. In the analysis of chargino 
exchange contributions it was found that the CP odd component of H2 is rather 
small while the analysis of H^ parallels the analysis of Hi with the only difference 
that the roles of the CP even and the CP odd components is reversed. Much the 
same situations occurs in this case and thus we omit the detailed discussion of 
these states. 



6 Conclusions 

In this paper we have developed a calculus for the derivatives of the eigen val- 
ues of the neutralino mass matrix with respect to the background fields which 
are in general dependent on CP violating phases. The calculus allows one to 
deduce the derivatives of the eigen values of the neutralino mass matrix analyti- 
cally even though the eigen values themselves cannot be gotten analytically in a 
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compact form. We use this calculus to obtain analytical results for the neutralino- 
Z-neutral Higgs exchange contribution to the masses and mixings in the CP even- 
CP odd neutral Higgs sector. The above computation along with the stop-top, the 
sbottom-bottom, the tau-stau and the chargino-W- charged Higgs exchange con- 
tribution computed previously provide us with a complete one loop contribution 
to the Higgs mass matrix with the inclusion of CP phases. This full one loop result 
was then used to discuss the phenomenon of CP violation in the neutral Higgs sec- 
tor. The numerical analysis shows that the mixings between the CP even and the 
CP odd sectors are significantly affected by the neutralino exchange contribution. 
The mixing of the CP even and the CP odd Higgs sector have many important 



consequences ||15|, [Tq , 0. Thus one consequence is that CP even-odd mixing af- 
fects the couplings of the Higgs bosons with quarks and leptons and this effect can 
be discerned in Higgs searches in collider experiments. Another important impli- 
cation is that the CP even-odd mixing will affect the relic density analysis and 
thus modify the parameter space allowed by the relic density constraints. Further, 
since the couplings of the quark and leptons with the Higgs are affected due to 
the CP even-odd mixing there will also be an effect of these mixings on detection 
rates in the direct searches for dark matter. It would be interesting to carry out 
an analysis of these phenomena. 
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7 Appendix A: Neutralino eigen values and deriva- 
tives 

The characteristic equation for the square of the neutralino mass is F(X)= DetiNlK 
M^o~XI) = where A represents the square of the neutralino mass eigen values. 
It can be expanded as 

F(A) = X^ + aX^ + bX^ + cX + d = (32) 



14 



In the above a,b,c and d are computed using Eq.(14). The computation of the 
co-efficients is done to leading and to next to the leading order in an expansion in 
M|/M| where Ms stands for the soft SUSY parameters. Thus, e.g., a is expanded 
to 0(M|) and 0(M|) orders (it is actually exact when expanded to this order), b 
is expanded to 0(M|) and 0{M'i;M'^) orders etc. The analysis for a, b and c (d 
does not enter in Eqs. (23-29) and is not exhibited) gives 

a = -[|mi|2 + |m2p + 2|/xp + 2M|] (33) 



U I ~ |2| ~ |2 I I |4 I ol |2/| ~ |2 I I ~ |2\ 

= \mi\ \m2\ + +2|/i| + |m2| j 

-FM|[|mi|^ + Imsl^ + 2|//|^ + (|mi|^ - |m2ncos26'vK 
— 4 cos /3 sin /3C^ 1 7712 1 1 A* I cos 72 — 4cos/3sin/35'^|mi||//| cos 71] (34) 

where = -rr^ and 5",^ = -^rhi 

^ 91+92 ^ 91+95 

ol |2| ~ |2| ~ |2 I |4/| ~ |2 I I ~ |2\ 

c——2\jj,\ \mi\ \m2\ — (|mi| + |m2| ) 
-F4M|sin/?cos/3|/x|[(|m2p + |//.n5'^|mi| cos 71 + {\mi\'^ + |/x|^)C^|m2| cos 72] 

(35) 

The derivatives dXi/d^a can be gotten exphcitly as follows: 

dXi _ a^A^ + b^X^ + CqA + da , 

4A3 + 3aA2 + 26A + c ^ ^ 

The second derivatives are given by 
d^X, ^ ^ Ja X^ + b X^ + c X + dX ^^^, ^ ^^^2 ^ + ^^^(,2A^ + 6aA + 26) 



d^ad^p ' (4A3 + 3aA2 + 26A + c) 

^ (a^A^^ + baX'^ + CaX + da) ^3^^_^2 ^ 26^A + c^) + ^^^'^^ ^ ^^^^ ^ ^ ^^"^ 



(4A3 + 3aA2 + 26A + c)2 ^ p ■ ■ ^4^3 + + 26A + c) 

(ttapX + 6q,/3A -f- Cq-aA -|- do 



2 



x(3a„A + 2baX + - (4;,3 + 3^^^ + 26A + c) ^'='^^^^^ 



where 



da d'^a 
9C ""^'^d^ 
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8 Appendix B: List of parameters 

The explicit evaluation of the co-efRcients ai,bi,ci, di is given below 

bi = -fl'2lA*ll"^2|t'2COS72 - tj'i|^||mi|f2COS7i 

+vi[\m,\'gl + Im^l'gl + {gl + gl)\fi\'] (39) 

ci = giilfhil'^ + |/ini/u||m2|v2Cos72 
+gl{\m2\'^ + \fi\^)\^i\\mi\v2 cos 71 
-g^\li\'^\mi\'^vi - gl\ii\'^\m2\'^vi (40) 



di = -5(2|/i|^|mip|m2|t;2 COS72 

-gl\lif\m2\^\mi\v2Cosji (41) 

The co-efRcients a2,b2,C2,d2 can be gotten from ai,6i,Ci,di with the following 
interchanges 

a2^ai{vi< ^^2), b2^bi{vi< > V2), C2 ^ Ci{vi < ^ ^2), d2 = di{vi < > V2){4:2) 

The co-efficients as, 63, C3, ds are given as follows 

^3 = -5'2|"7,2||A*|^^2sin72 - g'i|mi||/i|^;2sin7i (43) 



C3 = gi{\ml\'^ + \iJ,\'^)\m2\\n\v2smj2 + gl{\ml\^ + \ii\^)\mi\\ii\v2smji (44) 



ds = -g'2l"i.i|^|//|^|m2|'y2sin72 - gl\m2\'^\iJ,\^\rhi\v2sm^i (45) 
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The co-cfRcicnts 03,63, 03,^3 can be gotten from 03,63,03,^3 with the following 



interchanges 



4 = 03(^1 < >V2), 63 = 63(^^1 < ^^'2), C3 = 03(^1 < ^^2), d'^ = d3{vi< >V2){^Q) 

Aq and Bq are given by 

^0 = ^{gl + gi)v,{Ml - Mlo) cos2/3 + {gf + gl)v2{Ml + M%) sin2/3 (47) 

So = -2{gl + gl)v2{Ml - M%) cos 2(5 + {gl + y2')^;i(M| + M%) sin 2/3 (48) 



17 



Figure Captions 

Fig.l: Plot of the CP even component (f)i of Hi (upper curves) and the CP odd 
component ipio of Hi (lower curves) including the stop, sbottom, stau, chargino 
and neutralino sector contributions as a function of the scale Q. The common 
parameters are itla = 300, tan (3 = 15, mo = 100, nii = 500, = .4, ,^2 = -5, 
tto = -3, l^ol = 1- The curves with circles are for 6^ = 0.1 and with squares for 
6*^ = 0.2 where all masses are in GeV and all angles are in radians. 

Fig. 2: Plot of Ai3 including the stop, sbottom, stau, chargino and neutralino sec- 
tor contributions vs the U(l) gaugino phase ^1. The common input for all the 
curves are rriQ = 100, mi = 500, Ma = 300, \Ao\ = 1, ao = 0.3, ,^2 = 0.5 and 
Q = 320. The five curves correspond to the pairs of tan/? and 6^ values as follows. 
The curve with A13 = 301 at ,^1 = corresponds to tan/? = 5, 6^^ = .4. Similarly 
the curves with values of A13 = 406 at ^1 = correspond to tan/? = 6, = .6, 
Ai3 = 416 at ,^1 = correspond to tan/? = 10, 6^ = .2, A13 = 501 at ,^1 = cor- 
respond to tan/? = 8, = .8, and A13 = 579 at ,^1 = correspond to tan/? = 15, 
O/j, = .3 where all masses are in GeV and all angles are in radians. 

Fig. 3: Plot of A23 including the stop, sbottom, stau, chargino and neutralino sec- 
tor contributions vs the U(l) gaugino phase ^1 for the same input parameters as in 
Fig. 2. The curves with the same symbols as in Fig. 2 have the same common inputs. 

Fig. 4: Plot of the CP even component 0i of Hi (upper curves) and the CP odd 
component ipiD of Hi (lower curves) including the stop, sbottom, stau, chargino 
and neutralino sector contributions as a function of the U(l) gaugino phase ,^1 for 
the same inputs as in Fig.2. The curves with the same symbols as in Fig. 2 have 
the same common inputs. 

Fig. 5: Plot of the CP even component 0i of Hi (upper curves) and the CP odd 
component ipiD of Hi (lower curves) including the stop, sbottom, stau, chargino 
and neutralino sector contributions as a function of the ^2- The common parame- 
ters are: m^ = 300, Q = 320, mo = 100, mi = 500, ao = -3, |^o| = Ij = -4. For 
the curves with diamonds tan/? = 15, ^1 = 1.5, for squares tan/? = 8, .^1 = 1.5, for 
triangles tan/? = 8, ,^1 = 0.5, and for circles tan/? = 10, ^1 = 1.5 where all masses 
are in GeV and all angles are in radians. 



18 



Fig. 6: Plot of the CP even component (f)i of Hi (upper curves) and the CP odd 
component ipio of Hi (lower curves) including the stop, sbottom, stau, chargino 
and neutralino sector contributions as a function of 9^^. The common parameters 
are niA = 300, Q = 320, mo = 100, nii = 500, G = -5, = .3, \Ao\ = 1. For 

2 

curves with diamonds tan/9 = 15, = 1.5, for squares tan (3 = 8, C,i = 1-5, and 
for triangles tan/3 = 8, = 0.5 where all masses are in GeV and all angles are in 
radians. 

Fig. 7: Plot of the CP even component 0i of Hi (upper curves) and the CP odd 
component ipnj of Hi (lower curves) including the stop, sbottom, stau, chargino 
and neutralino sector contributions as a function of tan (3. The common input 
parameters for the curves are = 300, Q = 320, niQ = 100, mi = 500, = .5, 
^2 = -5, tto = -3, and \Ao\ = 1. For the curves with diamonds, 6^ = .4, for squares 
9fj, = .6, and for triangles 9^ = .8 where all masses are in GeV and all angles are 
in radians. 
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